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1. Introduction 

A fundamental problem in quantum computing and, more generally, in 
quantum information processing 1 is to experimentally achieve any theoret- 
ically designed quantum channel with a fixed device, being able to program 
the channel on the state of an ancilla. This problem is of relevance for 
example in proving the equivalence of cryptographic protocols, e. g. prov- 
ing the equivalence between a multi-round and a single-round quantum 
bit commitment 2 . What makes the problem of channel programmability 
non trivial is that exact universal programmability of channels is impossi- 
ble, as a consequence of a no-go theorem for programmability of unitary 
transformations by Nielsen and Chuang 3 . A similar situation occurs for 
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universal programmability of POVM's 4 ' 5 . It is still possible to achieve pro- 
grammability probabilistically 6 , or even deterministically 7 , though within 
some accuracy. Then, for the deterministic case, the problem is to deter- 
mine the most efficient programmability, namely the optimal dimension of 
the program-ancilla for given accuracy. Recently, it has been shown that 
a dimension increasing polynomially with precision is possible: however, 
even though this is a dramatical improvement compared to preliminary 
indications of an exponential grow 8 , still it is not optimal. 

In establishing the theoretical limits to state-programmability of chan- 
nels and POVM's the starting problem is to find the joint system-ancilla 
unitary which achieves the best accuracy for fixed dimension of the ancilla: 
this is exactly the problem that is addressed in the present paper. The 
problem turned out to be hard, even for low dimension, and here we will 
give a solution for the qubit case, for both system and ancilla. 

2. Statement of the problem 

We want to program the channel by a fixed device as follows 



with the system in the state p interacting with an ancilla in the state a via 
the unitary operator V of the programmable device (the state of the ancilla 
is the program). For fixed V the above map can be regarded as a linear map 
from the convex set of the ancilla states srf to the convex set of channels for 
the system ^ . We will denote by SPy^ the image of the ancilla states srf 
under such linear map: these are the programmable channels. According 
to the well known no-go theorem by Nielsen and Chuang it is impossible 
to program all unitary channels on the system with a single V and a finite- 
dimensional ancilla, namely the image convex C ^ is a proper subset 
of the whole convex c € of channels. This opens the following problem: 

Problem: For given dimension of the ancilla, find the unitary 
operators V that are the most efficient in programming channels, 
namely which minimize the largest distance e(V) of each channel 
C G c (> from the programmable set S?vM : 



As a definition of distance it would be most appropriate to use the CB- 
norm distance \\C — V\cb- However, this leads to a very hard problem. We 



VvAp) = ^2[V{p®a)V^] 



(1) 



e(V) = max min S(C, V) = max min S(C, Vv CT ). 



(2) 
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will use instead the following distance 

S(C,P) = y/l-F(C,P), (3) 

where FiC^V) denotes the Raginsky fidelity 9 , which for unitary map C = 
U = U ■ W is equivalent to the channel fidelity 1 

F(u,r) = ^\ Tr ^\^ ( 4 ) 

i 

where C = J2i Ci • C\ . Such fidelity is also related to the input-output 
fidelity averaged over all pure states F io (U, V), by the formula F io (U, V) = 
[1 + dF(U, V)]/ (d+ 1). Therefore, our optimal unitary V will maximize the 
fidelity 

F(V)= min F(U,V), F{U, V) = maxF(W, V v a) (5) 

(7GU(H) aesf 

3. Reducing the problem to an operator norm 

In the following we will use the GNS representation \^f)) = ® -01-0) of 
operators ^ e B(H), and denote by X 1 the transposed with respect to the 
cyclic vector |7)>, i. e. |*)) = (tf <g> 010) = i 1 ® * T )I0)> and b y x * thc 
complex conjugated operator X* = (X T )t, and write \v*) for the vector 
such that (|f)(f | ® 7)|0) = l w )l w *)- Upon spcctralizing thc unitary V as 
follows 

V = £e**|*fc»«*k|, (6) 
k 

we obtain the Kraus operators for the map Vv,a{p) 

TvAp) = J2 c -mpcL, c nm = J2^kK){v* m \^lV>^ (7) 

nm k 

where \v n ) denotes the eigenvector of a corresponding to the eigenvalue A„. 
We then obtain 

£ | Tr[Ct ro t/]| 2 = £ e^-"*) Tr[*t [/t^^T^t 

nm fc/i (8) 

= Tr[<7 T S(tf,V) t S(E7, V)] 

where 

s(t/,v) = ^ e - 4efc *^* fc . (9) 

k 

The fidelity (5) can then be rewritten as follows 

F(U,V) = ±\\S(U,V)f. (10) 
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4. Solution for the qubit case 

The operator S(U, V) in Eq. (9) can be written as follows 

S{U,V) =Tr 1 [(f/ T ®7)F*]. (11) 

Changing V by local unitary operators transforms S(U, V) in the following 
fashion 

S(U, (Wi <g> W 2 )V(W 3 <g> W 4 )) = W^W^t/W^, V)W 4 *, (12) 

namely the local unitaries do not change the minimum fidelity, since the 
unitaries on the ancilla just imply a different program state, whereas the 
unitaries on the system just imply that the minimum fidelity is achieved 
for a different unitary — say W\UW$ instead of U. 

For system and ancilla both two-dimensional, one can parameterize all 
possible joint unitary operators as follows 10 

V = (V^i(g)T^ 2 )exp[z(ai(Ti(g)CTi T +a2CT2«)CT2 T +a3CT3®CT3 T )](VF3(g)W4) • (13) 

A possible quantum circuit to achieve V in Eq. (13) can be designed using 
the identities 

[<r a ® (J a ,ap <g> up] = 0, 

C(a x <g> I)C = <7 x <Z>cr x , 

C(I®a z )C =-a z ®a z , ^ 
{e~^ az <g> e~T CT ^ C(a x <g> I)C (eT ,7 '®T'«j = Gy ® ayi 
where C denotes the controlled-NOT 

C=|0)<0|®/+|1)<1|®<t x . (15) 
This gives the quantum circuit in Fig. 1. The problem is now reduced to 
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Figure 1. Quantum circuit scheme for the general joint unitary operator V in Eq. (13). 
Here we use the notation = exp(i<j><JQ) with G = X,Y, Z. 



study only joint unitary operators of the form 

V = exp[(i(ai<7i <g> (Ti T + a 2 CT2 <8> ct 2 t + a 3 a 3 <g> cr 3 T )] . (16) 
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This has eigenvectors 



I*,-)) = -=\aj)), (17) 



1 

where aj, j — 0, 1, 2, 3 denote the Pauli matrices Co = I, &i = 02 = &y, 
cr 3 = u z . This means that we can rewrite S(U, V) in Eq. (9) as follows 

WV) = ^e-V";> (18) 

3=0 

with 

0o = ai+a 2 +a 3 , 9i = 2a» - 6q> • (19) 
The unitary U belongs to SU(2), and can be written in the Bloch form 

U = n I + in-<r, (20) 
with rik € R and n 2 , + |n| 2 = 1. Using the identity 

ajcrio-j = ejiai, e j0 = e 3j = 1, = -1,1 ^ 0, j, (21) 
we can rewrite 

S(U, V) = hoi + n ■ er, (22) 

where 



3 



hj =tjrij , < j < 3, to = — e 



(23) 



(24) 



3=0 

tj =er l0a + e~ i9i -t , l<j< 3, tj = \tj\e 1 ^ , < j < 3, 

It is now easy to evaluate the operator S(U, V)^ S(U, V). One has 

S(U,V?S(U,V) =v I + v-<r, 

v a = \h Q \ 2 + \n\ 2 , v = i [2Q(n n*) + h* x n] . 

Now, the maximum eigenvalue of S(U, V)^ S(U, V) is vo + \v\, and one has 
\v\ 2 = E N 2 |n,| 2 -nfn| = 2^ |n,| 2 |^| 2 sin 2 (^ - (25) 

i,j=0 i,j=0 

whence the norm of S"(?7, V) is given by 

3 



\S{U,V)f =Y,n%\ 2 + 

3=0 



3 

2 £ n 2 n 2 |t 4 |%| 2 sin 2 (<^ _ . (26) 

ij=0 
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Notice that the unitary U which is programmed with minimum fidelity in 
general will not not be unique, since the expression for the fidelity depends 
on {n 2 }. Notice also that using the decomposition in Eq. (13) the minimum 
fidelity just depends on the phases {9j}, and the local unitaries will appear 
only in the definitions of the optimal program state and of the worstly 
approximated unitary. It is convenient to write Eq. (26) as follows 

ut + Vu-Tu. (27) 

= (M 2 , 

— (pj). One has the bounds 
u ■ t + \Fu 



where u 
\U\ 2 \t, ' 



\\S(U,V)f 
{nl,nl,nl,n%), t 



|ti| 2 ,M 2 ,N 2 ), and Ta = 



■ sin 



Tu > u ■ t > min \tj 



(28) 



and the bound is achieved on one of the for extremal points ui — Sij of 
the domain of u which is the convex set {it, Uj > 0, J2j u j = 1} (t ne 
positive octant of the unit four dimensional ball S+). Therefore, the fidelity 
minimized over all unitaries is given by 



nv) = ^vf\tj\ 



(29) 



The optimal unitary V is now obtained by maximizing F(V). We need 
then to consider the decomposition Eq. (13), and then to maximize the 
minimum among the four eigenvalues of S(U, V^SiU, V). Notice that tj = 
V ffj (1 e l8 '' I where H is the Hadamard matrix 



H = 



/I 
1 
1 

M 



i 
i 

-l 
-l 



l 



i/ 



(30) 



which is unitary, and consequently ^ . |t 3 | 2 = J2j |e 49j | 2 = 4. This implies 
that min 3 \tj\ < 1. We now provide a choice of phases 6j such that \tj\ = 1 
for all j, achieving the maximum fidelity allowed. For instance, we can 
take O =0,01 = ir/2,02 = n, 6% = ir/2, corresponding to the eigenvalues 
i, 1, —i, 1 for V. Another solution is O = 0, 0\ = — ir/2, 9 2 = ir, 63 = —tt/2. 
Also one can set Q{ — » — 9i. The eigenvalues of S(U, VyS(U, V) are then 
1,1,1,1, while for the fidelity we have 

nv) = ± = \, (3D 



max 
yeu(H® 2 ) 



and the corresponding optimal V has the form 



V = exp 



7T . 



•a x ± a z (g) a z ) 



(32) 
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A possible circuit scheme for the optimal V is given in Fig. 2. 



X ±1 



-e- 



Figure 2. Quantum circuit scheme for the optimal unitary operator V in Eq. (31). For 
the notation see Fig. 1. For the derivation of the circuit see Eqs. (14). 



We now show that such fidelity cannot be achieved by any V of the 
controllcd-unitary form 

2 

V = J2 V k®\iPk)(il>k\, (V'i|^2)=0, V U V 2 unitary on H ~ C 2 . (33) 
fc=i 

For spectral decomposition Vk = Y^j^i^ 9 *^^^)^^ the eigenvectors 
of V are \^jkfj = |<?^)|V'fc); an d tne corresponding operators are ^jk = 
\<j>j k ^)(ipl\, namely the operator S(U, V) is 

S(U,V) = ^e- ie i j) m)(^ k) \M k) )(r k \ , (34) 

with singular values £ 2 =1 e~ ie ™ {<f>f ] \U\<i)f } ) = TT [Vk U ]- Thcn ' thc op- 
timal program state is \iph), with h = argmax^ |Tr[Vj[/]|, and the corre- 
sponding fidelity is 

F{U,V) = \\Tr[V^U]\\ (35) 

and one has 

F(V) = mm F(U, V) = 0, (36) 

since for any couple of unitaries Vk there always exists a unitary U such that 
Tr[V^ £7] = for k = 1,2. Indeed, writing the unitaries in the Bloch form 
(20), their Hilbcrt- Schmidt scalar is equal to the euclidean scalar product 
in R 4 of their corresponding vectors, whence it is always possible to find a 
vector orthogonal to any given couple in M 4 . The corresponding U is then 
orthogonal to both Vk , and the minimum fidelity for any controlled-unitary 
is zero. 
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